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1. Motivation and objectives 

Study of turbulent flows in rotating reference frames has long been an area of 
considerable scientific and engineering interest. Because of its importance, the sub- 
ject of turbulence in rotating reference frames has motivated over the years a large 
number of theoretical, experimental, and computational studies (e.g., Greenspan 
1968, Bardina et al 1985, Jacquin et al 1990, Mansour et al 1991). The bulk of 
these previous works has served to demonstrate that the effect of system rotation 
on turbulence is subtle and remains exceedingly difficult to predict. 

A rotating flow of particular interest in many studies, including the present work, 
is examination of the effect of solid-body rotation on an initially isotropic turbulent 
flow. One of the principal reasons for the interest in this flow is that it represents the 
most basic turbulent flow whose structure is altered by system rotation but without 
the complicating effects introduced by mean strains or flow inhomogeneities. The 
assumption of statistical homogeneity considerably simplifies analysis and compu- 
tation. 

For an initially isotropic turbulence, it is well known that system rotation inhibits 
the non-linear cascade of energy from large to small scales. This effect is manifest in 
a reduction of the turbulence dissipation rate and associated decrease in the decay 
rate of turbulence kinetic energy (e.g., see Traugott 1958, Veeravalli 1991, Mansour 
et al 1992). An issue considerably less resolved, however, is the development of 
a two-dimensional state in rotating homogeneous turbulence. Both computations 
and experiments have noted an increase in integral length scales along the rotation 
axis relative to those in non-rotating turbulence (Bardina et al 1985, Jacquin et 
al 1990). Increase in the integral length scales has been thought to be a prelude 
to a Taylor-Proudman reorganization to two-dimensional turbulence. However, it 
has also been shown using direct numerical simulation (DNS) of rotating isotropic 
turbulence (Speziale et al 1987, Mansour et al 1992) that in the limit of very rapid 
rotation, turbulence remains isotropic and three-dimensional. In fact, Mansour et 
al (1992) showed that the evolution of rapidly rotating turbulence was accurately 
predicted using rapid distortion theory (RDT). Furthermore, Mansour et al also 
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showed that the RDT solution violates a necessary condition for occurrence of a 
Taylor-Proudman reorganization. 

It is worth noting that as is typically the case with DNS, the computations 
performed by Speziale et al. (1987) and Mansour et al. (1991,1992) were performed 
at low Reynolds numbers. In rotating turbulence at low Reynolds number, the 
effects of viscous decay progressively reduce the Rossby number and drive the flow to 
the RDT limit. Thus, other mechanisms for obtaining two-dimensional turbulence, 
e.g., through non-linear interactions which occur on a turbulence time scale, are 
precluded using DNS, and evolution of a two-dimensional state, therefore, requires 
significantly higher Reynolds numbers than can be attained using DNS. 

An issue closely connected to development of two-dimensional turbulence in rotat- 
ing flows is the existence of asymptotic self-similar states. The issue of self-similarity 
is a topic central to studies of turbulent flows (e.g., see Chasnov 1993). A similar- 
ity state is characterized by the predictability of future flow statistics from current 
values by a simple rescaling of the statistics; the rescaling is typically based on 
a dimensional invariant of the flow. Knowledge of the existence of an asymptotic 
similarity state allows a prediction of the ultimate statistical evolution of the flow 
without detailed knowledge of the complex, and not well understood, non-linear 
transfer processes. 

Large-eddy simulation (LES) is ideally suited for examination of the long-time 
evolution of rotating turbulence since it circumvents the Reynolds number restric- 
tion of DNS. The drawback is, of course, that it requires use of a model to parame- 
terize subgrid-scale stresses. However, large-scale statistics are relatively insensitive 
to the exact form of the model, and alternative approaches, i.e., laboratory experi- 
ments or direct simulations, are simply not feasible for examination of the long-time 
evolution of rotating flows. 

The principal objective of the present study has thus been to examine the asymp- 
totic state of solid-body rotation applied to an initially isotropic, high Reynolds 
number turbulent flow. Of particular interest has been to determine (1) the degree of 
two-dimensionalization and (2) the existence of asymptotic self-similar states in ho- 
mogeneous rotating turbulence. As shown in §2, development of a two-dimensional 
state is very pronounced; much more so than observed in previous studies using 
DNS. It is also shown that long-time evolution of quantities such as turbulence ki- 
netic energy and integral length scales are accurately predicted using simple scaling 
arguments. 

2. Accomplishments 

2.1 Simulation overview 

In the present study, the filtered Navier-Stokes equations for an incompressible 
fluid were solved in a rotating reference frame: 

V • u = 0 (1) 


d 1 

-f u • Vu = — Vp + i/ e V 2 u — 2Hxu. 

at p 


(2) 
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In (1) and (2), u is the velocity vector, p and p the fluid pressure and density, 
respectively, and ft is the rotation vector. For purposes of discussion, the rotation 
vector is considered to act along the z or “vertical” axis, fl = (0,0, ft). An eddy 
viscosity hypothesis was used to parameterize the subgrid- scale stresses. In this 
work, the spectral eddy viscosity of Chollet & Lesieur (1981) was modified for 
rotating turbulence 

V t = I'eo /(a) (3) 

where i/ e o is the “baseline” viscosity and /(a) a function accounting for the reduction 
of v e by system rotation. The baseline viscosity, v e o, from Chollet Lesieur, is 


i/ e 0 (A:|fc m ,t) 


0.145 + 5.01exp 


/ — 3.03fc m 

V * 



E(k m ,t ) 



(4) 


where k m is the maximum wavenumber magnitude of the simulation and E(k, t) 
is the spherically integrated three-dimensional Fourier transform of the co- variance 
|(u,(x, t)u;(x + r, t)) ((•) denotes an ensemble or volume average). The reduction 
in i/ e is expressed using f(a) 


f(a) = 


2 [(1 + a 2 ) 3 / 2 - a 3 - 1] 


where 


3a 2 

8ft 2 


a = 


3 E{k m )k* m 

(Cambon, private communication). 

The initial energy spectrum of the simulations was of the form 


[-Hi)'] 


where s is equal to 2 or 4, C 3 is given by 


(5) 

(6) 


(7) 


C 9 = 


j~2 s$ (3 + l) 

V 7T 1 * 3 ‘ * (s — 1) 


( 8 ) 


and k p is the wavenumber at which the initial energy spectrum is maximum. In 
this study, simulations with s ~2 and s = 4 were performed, corresponding to the 
initial energy spectra with a low wavenumber form proportional to either k 2 or fc 4 . 

Because the principal interest of this work was examination of the long-time evo- 
lution of rotating turbulence, it was necessary to use as large a value of k p as possible 
in order that flow evolution not be adversely affected by the periodic boundary con- 
ditions used in the simulations; adverse affects occurring when the integral length 
scales of the flow become comparable to the box size. Another important con- 
sideration in these simulations was the aspect ratio of the computational domain. 
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Because of the rapid growth of turbulence length scales along the direction of the 
rotation axis, it was necessary to use a computational box which was longer along 
the rotation axis than in the other directions. Preliminary calculations of rotating 
turbulence on cubic domains had shown an adverse affect of periodicity because 
of the rapid integral scale growth along the rotation axis. Numerical experiments 
showed that it was necessary to use to a computational box which was four times 
larger along the rotation axis than in the directions orthogonal to the rotation vec- 
tor. Four times as many collocation points were used in the vertical direction in 
order to avoid any effects of grid anisotropy at the smallest resolved scales. 

Simulations were performed using resolutions of 96 x 96 x 384 and 128x128x512 
collocation points. The governing equations (1) and (2) were solved using the 
pseudo-spectral method developed by Rogallo (1981). The statistical evolution of 
the flow using either resolution was essentially the same, and, therefore, only the 
results from the higher resolution computation are reported in this brief. The 
maximum physical wavenumber of the 128 x 128 x 512 computations was 95 (for a 
computational domain having total volume of 87r 3 ). The initial root-mean-square 
velocity fluctuation uq in (7) was equal to unity and k p = 75. Following Chasnov 
(1993), the initial energy spectrum was set to zero for wavenumbers greater than 
93 to allow the sub grid-scale eddy viscosity to build up from zero values. For each 
spectrum type, i.e., low wavenumber part proportional to k 2 or fc 4 , simulations were 
performed with Q = 0, 0.5, and 1.0. 

2.2 Results 

The instantaneous power-law exponent (i.e., the logarithmic time derivative) of 
the mean-square velocity fluctuation, (u 2 ), is shown in Figure 1 for each rotation 
rate and spectrum type used in the simulations. The time axis in Figure 1 and 
following figures has been made dimensionless using the eddy turnover time in the 
initial field 

r(0) = L u (Q)/(u 2 )'' 2 (9) 

where T u (t) is the velocity integral scale defined at time t as 


r ^ 7T f 0 °°k-'E(k,t)dk 
“ () ~2 / 0 °° E(k,t)dk ' 

In isotropic turbulence, L u is two-thirds the usual longitudinal integral scale mea- 
sured in experiments (see also Chasnov 1993). Throughout this work, u k 2 spectrum” 
refers to an initial energy spectrum E(k) with low wavenumber part proportional 
to k 2 while “A: 4 spectrum” refers to an initial E(k) with low wavenumbers propor- 
tional to fc 4 . For both spectrum types, the characteristic effect of system rotation 
in reducing the decay rate of (u 2 ) is evident in Figure 1. Furthermore, it is also 
clear from the Figure that the value of the decay exponent for non-zero rotation 
rates at long times is independent of fi, depending only on the form of the initial 
energy spectrum. 

It is possible to predict the exponents in Figure 1 if one assumes that the asymp- 
totic scaling of (u 2 ) is dependent on the form of E(k) at low wavenumbers and 
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FIGURE 1. Time development of the power law exponent of (u 2 ) in rotating 

turbulence. ©■ ■ o ,11=0 , Ar 2 spectrum; , ft = 0.5, Ar 2 spectrum; , 

ft = 1.0, k 2 spectrum; , ft = 0, k 4 spectrum; , ft = 0.5, k 4 spectrum; 

, ft = 1.0, Ar 4 spectrum. 

independent of viscosity. For high Reynolds number turbulence, this is a reasonable 
assumption since the direct effect of viscosity occurs at much larger wavenumber 
magnitudes than those scales which contain most of the energy. Thus, it is possible 
to derive expressions for the asymptotic scaling of (u 2 ) using an expansion of the 
energy spectra near k = 0: 

E(k) = 2irk 2 (A Q +A 2 k 2 + ...) (11) 

where Ao. A 2 , ... are the Taylor series coefficients of the expansion. As shown 
by Batchelor & Proudman (1956), assuming that the velocity correlation tensor 
(u,(x)tij(x + r)) is analytic at k = 0 results in Ao = 0 and a time-dependent value 
of A 2 . On the other hand, SafFman (1967) showed that it is also physically possible 
to create an isotropic turbulence with a non- zero value of Ao, which is also invariant 
in time. As also shown by Chasnov (1993) the asymptotic scalings of (u 2 ) for these 
two cases are 

(u 2 ) oc A^ 5 f _6//5 (Ar 2 spectrum) (12) 

and 

(u 2 ) a A^ 7 * -10 / 7 (Ar 4 spectrum). (13) 

It may be observed from Figure 1 that the agreement between the LES results and 
(12) and (13) is excellent for both spectrum types. An analysis similar to that 
leading to (12) and (13) may also be performed for flows having non-zero rotation 
rates. For non-zero ft, the asymptotic scalings of (u 2 ) are predicted to be 

(u 2 ) a A^t 3 / 5 ft 3 / 5 (Ar 2 spectrum) 


(14) 
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(u 2 > oc Al /7 t 5 / 7 fl 5 / 7 (A: 4 spectrum). (15) 

For the k 2 spectrum, the value of the decay exponent from the rotating flows in 
the asymptotic region, approximately -0.64, is in very good agreement with (14). 
Similar agreement between the predicted value of the decay exponent, -5/7, and the 
measured values for the A: 4 spectrum is also observed. 

Evolution of the integral length scales are shown in Figures 2a and 2b. The 
vertical integral scale is defined as 


L v - J (ui{x,y,z)ui(x,y,z + r))dr 


while the horizontal integral scale is given by 


L h 


1 [ (uj(x)Mx + r)) 

(u 2 ) J 2nr 


(16) 


(17) 


Figure 2 clearly shows the significantly greater growth in time of the vertical integral 
scales relative to their horizontal counterparts in rotating turbulence. Also shown 
in the Figure is the velocity integral scale, L u , for ft = 0. It is evident from Figure 2 
that the horizontal integral scales in the simulations with non-zero ft are essentially 
independent of rotation rate and evolve similarly to the length scale from the non- 
rotating case. 

The results in Figure 2 may be used to deduce a posteriori the asymptotic scaling 
laws of the integral scales. For the A: 2 spectrum, dimensional arguments and the 
LES results in Figure 2a give the following dependence of the length scales on the 
invariant Ao, and ft 


Lh oc A l J 5 t 2 / 5 (fc 2 spectrum), (18) 

i.e., no dependence of Lh on fh The appropriate scaling of the vertical length scales 
for the k 2 spectrum is 


L v oc (Ar 2 spectrum) (19) 

since the long-time growth of L v is observed from the LES results to be directly 
proportional to time. Similarly, dimensional arguments together with the results in 
Figure 2b can be used to deduce the length scale dependence on A 21 t, and f2 for 
the k 4 spectrum: 

Lh oc Al^t 2 / 7 ( A* 4 spectrum), (20) 

similar to the non-rotating case. For the vertical length scales 

L v oc A l 2 ^ 7 tti 5 / 7 ( k 4 spectrum). 


( 21 ) 
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FIGURE 2. Time development of the integral length scales in rotating turbulence. 

(a) k 2 spectrum, ( b ) k 4 spectrum. , L u , ft = 0; , Z*, ft = 0.5; , 

L v , ft = 0.5; , Efiy 0 = 1.0; , 17 — 1.0. 


As was the case for the simulations possessing a k 2 spectrum, long-time evolution 
of L v in the simulations having a k 4 spectrum was also observed to be directly 
proportional to time. 

As shown in Figures 2a and 2b, the evolution of the flow in the direction along 
the vertical axis is strongly enhanced relative to the horizontal directions. Rapid 
growth of the vertical length scales provides an indication of evolution towards a 
two-dimensional state. This can be more clearly seen through examination of the 
energy spectrum as a function of spatial wavenumber k as well as the cosine of the 
polar angle in wave space (schematically illustrated in Figure 3). 

Shown in Figure 4a is the energy spectrum as a function of both k and cos 6 from a 
simulation with ft = 0 and possessing an initial spectrum with low wavenumber part 
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FIGURE 3. Wavenumber space showing rotation vector and polar angle #. For 
simulations of rotating turbulence, the pole is defined as # = 0 (cos# = 1); the 
equator as 9 = n/2 (cos# = 0). Energy and transfer spectra are considered to be 
functions of both k and cos #. 

proportional to fc 4 . It may be observed from the Figure that, as expected, the energy 
is essentially equi-partitioned with respect to cos#. Plotted in Figure 4b is the 
energy spectrum from a simulation with ft = 1 (and k 4 spectrum). It is clear there 
is a marked concentration of energy in the equatorial plane, # = 7t/2; the Figure 
provides very strong evidence of the development of two-dimensional turbulence. 
This result is also in sharp contrast to the previous examinations of E(k,co$6) by 
Mansour et al. (1992) using direct numerical simulation. Mansour et al. found only 
a slight tendency for a concentration of energy near the equator. Because of viscous 
dissipation in their simulations, it was not possible for Mansour et al to integrate the 
flow fields for long enough times in order to observe development of two-dimensional 
turbulence. It is important to emphasize that development of a two-dimensional 
state as demonstrated by Figure 4b cannot be captured by DNS because of viscous 
decay. LES circumvents this restriction and permits long enough integrations such 
that the non-linear interactions responsible for two-dimensionalization can occur. 

Further evidence of the profound effect of rotation is contained in Figures 5a 
and 5b. Figure 5a is the transfer function, T(k, cos#), from a simulation with 
k 2 spectrum and Q = 0. The Figure shows the expected behavior, i.e., negative 
transfer at low wavenumbers and positive transfer at higher wavenumbers. It is also 
reasonably clear from Figure 5a that, as was the case with the energy spectrum, the 
transfer term is independent of #. The transfer term from the simulation with a k 2 
spectrum and fl = 1.0 is shown in Figure 5b. As is clear from the Figure, rotation 
has substantially altered the transfer spectrum. For values of cos# near 1 (the pole 
in i-space), the energy transfer is small for the low wavenumbers and zero at the 
higher wavenumbers. For cos# near 0 (the equator in fc-space), the transfer term 
is actually positive at low wavenumbers, indicating a transfer of energy into these 
modes. 
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k 


FIGURE 4. Energy spectrum as a function of wavenumber k and cosine of the polar 
angle 6; spectrum obtained from LES with k 4 spectrum, (a) fl = 0 at t/r( 0) = 427, 
( b ) Q = 1.0 at t/r( 0) = 575. 


Time development of the two- and three-dimensional components of the kinetic 
energy are shown in Figures 6a and 6b for both spectrum types and each rotation 
rate. The two-dimensional component of the energy is obtained from Fourier modes 
in the plane k z = 0 while the three-dimensional component is from Fourier modes 
with k z / 0. The behavior is similar in both Figures and corroborates many of 
the aspects of the flows observed in the previous Figures. As expected, it may be 
observed that the decay of the energy is reduced with increasing rotation rates. 
More importantly, the Figures also show that for non- zero £2 the two-dimensional 
energy actually increases at later times in the flow evolution, consistent with Figure 
5b showing a transfer of energy into the low wavenumber modes in the equatorial 
plane. 

Figures 7a and 7b show the temporal evolution of the diagonal components of 
the anisotropy tensor of the Reynolds stress 


_ <«««>) &ij 

tJ (u 2 ) 3 


( 22 ) 


for each spectrum type and £2 == 0 and 1.0. It is interesting to note that, while 
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FIGURE 5. Transfer spectrum as a function of wavenumber k and cosine of 
the polar angle 6 ; spectrum obtained from LES with k 2 spectrum, (a) Q = 0 at 
i/r( 0) = 725, (6) Q = 1.0 at t/r( 0) = 552. 


the flow is becoming two-dimensional under the influence of system rotation, the 
distribution of kinetic energy amongst the three components shows little departure 
from isotropic values. As shown in the Figures, the vertical fluctuations are slightly 
enhanced by rotation relative to the fluctuations in the horizontal plane. The de- 
velopment of the diagonal components of the anisotropy tensor of the vorticity 



is shown in Figure 8. In contrast to the Reynolds stress anisotropy, this Figure 
shows an enhancement of vertical vorticity relative to the horizontal components. 

It is worthwhile to point out that recent work by Bartello et al (1993) showed 
a much stronger departure from isotropy of the component energies in rotating 
turbulence. Bartello et al. used cubic domains to examine the long-time evolution 
of the flow and consequently the evolution of the component energies was adversely 
impacted. Based on their findings, Bartello et al. concluded that the long-time 
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FIGURE 6. Evolution of two-dimensional and three-dimensional energy in LES. 

(a) k 2 spectrum, ( b ) k* spectrum. , 2D energy, ft = 0; o o , 3D energy, 

= 0; , 2D energy, ft = 0.5; , 3D energy, ft = 0.5; , 2D energy, 

Q = 1.0; , 3D energy, Q — 1.0. 


state of homogeneous rotating turbulence was two-dimensional and two-component. 
However, it is clear from Figures 7 and 8 that the asymptotic state of homogeneous 
rotating turbulence is two-dimensional but three- component. 

3. Future plans 

LES results presented in this work demonstrated the existence of asymptotically 
self-similar states in rotating homogeneous turbulence. Additional investigations 
are planned to further corroborate this finding, e.g., examining the collapse of the 
spatial spectra under the appropriate scalings. Related to this issue are the par- 
ticular scalings found in the present work. For example, the scaling laws for (u^) 
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t/r(Q) 

FIGURE 7. Evolution of Reynolds stress anisotropy tensor, (a) k 2 spectrum, (b) 
k 4 spectrum. Note that the curves for ft = 0 are essentially coincident and lie along 
the axis b aa — 0. , i>n, ft = 0; o - -e , b^ 2 , = 0; , f>33! = 0; , 

b lu Q = 1.0; , 622, ft = 10; , 633, ft = 10. 

in rotating turbulence (Equations 14 and 15) can be obtained using the transport 
equation for (u 2 ) together with dimensional analysis. The scaling laws for the 
length scales, however, (Equations 18-21) rely on only dimensional analysis and a 
posteriori examination of LES results. Development of a more rigorous approach 
for predicting the length scale evolution observed in the present work would be 
desirable. 

This study also showed the utility of large-eddy simulation. As shown in this 
work, development of an asymptotic state which is two-dimensional (but three- 
component) is very significant, much more so than in previous examinations of 
rotating turbulence using DNS. Given the encouraging results obtained to date, a 
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FIGURE 8. Evolution of vorticity anisotropy tensor, (a) k 2 spectrum, (6) Jb 4 
spectrum. Note that the curves for = 0 are essentially coincident and lie along 

the axis v aQ[ — 0 . , ^ = 0 ; o © , V221 = 0 j , V33, H = Oj , 

t'n, f) = 1.0; , ^22, ^ = 1.0; , U33, ft = 1.0. 


likely extension of this work will be use of the LES database for developing a high 
Reynolds number extension of the k-e model for rotating flows. 
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